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Gravitational perturbations about a Kerr black hole in the Newman-Penrose formalism are 
concisely described by the Teukolsky equation. New numerical methods for studying the evolution 
of such perturbations require not only the construction of appropriate initial data to describe the 
collision of two orbiting black holes, but also to know how such new data must be imposed into the 
, Teukolsky equation. In this paper we show how Cauchy data can be incorporated explicitly into 

the Teukolsky equation for non-rotating black holes. The Teukolsky function ^ and its first time 
derivative dt'^ can be written in terms of only the 3-geometry and the extrinsic curvature in a gauge 
invariant way. Taking a Laplace transform of the Teukolsky equation incorporates initial data as a 
source term. We show that for astrophysical data the straightforward Green function method leads 
to divergent integrals that can be regularized like for the case of a source generated by a particle 
coming from infinity. 
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■ I. INTRODUCTION AND OVERVIEW 

oo 

' Black holes coalescence is considered to be one of the strongest astrophysical sources of gravitational radiation and a 
I , primary candidate to be detected by gravitational wave observatories now under construction. For this reason, several 
' groups of researchers are now attacking the problem of solving Einstein's equations numerically, with supercomputers 
[ — _ [Q. Many difficulties must still be solved in this approach, such as the presence of instabilities in the numerical 
On ■ evolution due to the nonlinearity of the Einstein's equations and that of finding a better qualitative prescription for 
' truly astrophysical initial data representing black holes collisions. 

Meanwhile, the perturbation theory represents a valid and complementary approach to full numerical simulations, 
' with some clear advantages with respect the supercomputer project: It is a considerably more economical approach and 
(2J[), it is semianalytical. This latter aspect is a very important one since it allows to study some fundamental and conceptual 
1!V ■ problems. Among the new theoretical results in perturbation theory stands out the close limit approximation which 
. ^ \ approximates the collision of two black holes by a single perturbed one. The subsequent extension to second order 
^ • perturbations , and to moving (head-on) holes Q , confirmed the success of this approximation. Also recently [^-Q , 
5—1 ' the collision of a small and a big black hole by perturbative methods have been revisited in order to incorporate non- 
vanishing initial data. However, only perturbations about Schwarzschild black holes have been studied so far, by use 
of the Zerilli-Moncrief equation (a gauge invariant description). The success of these techniques encourages now 
to extend them to the more realistic rotating background. In particular, for the most plausible astrophysical scenario 
of a final single rotating hole, the first order "curvature" perturbations are compactly described by the Teukolsky 



equation 1 10 ll|] , which in the non-rotating case reduces to the so called Bardeen-Press equation ]lz 

The Teukolsky equation has been studied already in the early seventies in the frequency domain. In order to avoid 
the important problem of the imposition of the initial data into that equation, the computation of the gravitational 
radiation have been carried out in the case of unbounded particle trajectories (or circular motion) p3[ |. The divergent 
integrals encountered when one wants to calculate the gravitational radiation generated by a test mass falling into a 
black hole from infinity using the standard Green function, have been made finite by systematically discarding infinite 
surface terms [Q. In a different way to tackle this problem, Nakamura and Sasaki [|l^ transformed, in the frequency 
domain, the Teukolsky equation into other "more regular" equation. Only very recently, Poisson |]l6| (in the non- 
rotating case) and Campanelli and Lousto (for the rotating hole) showed that there is nothing intrinsically wrong 
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with the radial Teukolsky equation when deahng with unbounded source terms, and that the divergent integrals can 
be regularized in a natural way. In the Appendix, we will show how those results about regularization also apply 
when dealing with initial data for a bounded problem (like for the close limit case) . 

The problem of imposition of initial data in the Teukolsky equation is not, of course, restricted to the frequency 
domain form of the equation. Very recently an evolution code to integrate the Teukolsky equation in its time-domain 
form have been developed fl^ . However, until now the evolution of initial data have been restricted to the simple 
case of a bell-shaped burst, which cannot represent accurately realistic astrophysical initial data for the late stages of 
binary black hole coalescence. The imposition of the initial data into the Teukolsky equation, to our knowledge, seems 
not to have a clear framework yet. Consequently, in this paper we address to the important issue of the imposition 
of the initial data into the Teukolsky equation. In the non-rotating case, the gauge invariant Teukolsky functions ■00 
(ingoing) and -04 (outgoing), and their first time derivatives can be completely expressed in terms of the Moncrief 
waveform (and its time derivative), and then, in terms of the 3-geometry and the extrinsic curvature. By evaluating 
these functions on a constant time hypersurface for any specific astrophysical initial data, one is in principle able to 
numerically evolve the Teukolsky equation and to study very interesting astrophysical scenarios, like the close limit 
approximation [jl^ . The imposition of initial data into the Bardeen-Press equation can suggest what features 
need to be generalized in the rotating case. 

The paper is organized as follows: Section II is devoted to the problem of expressing the Teukolsky functions ipo 
and -04 on a given hypersurface in terms of the 3-geometry and the extrinsic curvature of that hypersurface. The 
computations are specifically performed in the Regge- Wheeler gauge j20|] for both even and odd parity waves and then 
reexpressed in a gauge invariant way by use of the Moncrief functions. In Sec. Ill we discuss the applicability of these 
relations and the possibility of generalization of these results to the rotating case. In the Appendix we briefly review 
the Teukolsky equation in its frequency domain and use the Laplace transform to mathematically incorporate initial 
data into the radial equation as an effective source term. We explicitly replace Brill-Lindquist initial data into the 
source term of the radial Teukolsky equation. The source term is shown to lead to a divergent integral as the domain 
of integration in the coordinate r extends to infinity. The regularization procedure already used for unbounded sources 
is employed successfully in the present case as described in the Appendix. 

The notation and conventions employed in this paper are those of Misner, Thorne and Wheeler ||2^. That is, 
the metric signature is (— -f + +); the Riemann tensor is defined by the equation R^f^^s — ^"/s-y s ~ ^"iss j + 
F'^^^F"^^ — V^pgY"^^ and the Ricci tensor by equation i?^^ = R^^^u- Throughout the paper an overdot indicates 
differentiation with respect to the time variable t, a prime differentiation with respect to the radial coordinate r, an 
overbar represents complex conjugation, and units are such that G = c = 1. The conventions for the NP quantities, 
i. e. spin coefficients, operators (denoted by an overhat and the Kinnersley tetrad), employed in this paper are those 
given in the table of Ref . [|l^ . 



II. INITIAL DATA FOR THE TEUKOLSKY EQUATION 



In the Newman-Penrose formalism, gravitational perturbations of the Kerr metric can be decoupled and described 
by the Teukolsky [0| equation. In the Boyer-Lindquist coordinates (t, r, 9, ip) this equation takes the following form 
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(2.1) 



where M is the mass of the black hole, a its angular momentum per unit mass, E = + cos^ 9, A = — 2M r + a^, 
and s is the spin parameter characterizing the perturbation. 
The field (where p — l/{r — ia cos 9) ) 



'^{t,r,9,(p) 



p ^^J^l'' = -p ^Cnrnnrn for S = -2 



for s = +2 



(2.2) 



represents either the outgoing radiative part of the perturbed Weyl tensor, ijj'^^ (s 
^„ _ , T-i,„ a;„„ „/,(") „/,(") 



-2), or the ingoing radiative 



part, {s = +2). The corresponding unperturbed quantities and -^q vanish identically in the Kerr spacetime. 
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The knowledge of only one of these two first order complex NP components, formed by projecting the perturbed Weyl 
tensor along the Kinnersley tetrad, is sufficient to uniquely and completely specify the gravitational perturbation 
(up to changes of the Kerr parameters a and M). In fact, ip'i'^ and ipg^"^ are gauge invariant under both coordinate 
transformations and infinitesimal tetrad (Lorentz) rotations. These two scalars carry information, in their real and 
imaginary parts, about the two dynamical degrees of freedom of the perturbed field. More precisely, in the nonrotating 
case, the real part of \1/ contains information about the even parity (or polar) perturbations of the metric (which are 
invariant under the transformation ip ^ —tp) while its imaginary part contains information of odd parity (or axial) 
perturbations of the metric (that change sign under a sign reversal of ip) . The source term T in Eq. ( |2.lD is constructed 
out of the Tl}i^ and its expression can be found in Refs. [ pl] , p7[ . 

In order to study perturbations around a Kerr hole one has to supplement the Teukolsky equation with initial 
data that allow to start the integration. So far the Teukolsky equation has been studied (in the frequency domain) 
for perturbations generated by particles coming from infinity, i.e. with vanishing initial data. There are, though, 
astrophysically relevant situations where initial data contribution are important. This is the case of the very successful 
close limit approximation We can divide the initial value problem for the Teukolsky equation into two subproblems: 
First, one must find what Cauchy data can better represent astrophysical scenarios for colliding black holes. Slices of 
constant Boyer-Linquist time of the Kerr geometry give a non-conformally flat three-geometry, making it incompatible 
with all initial data studied so far. It is not known wheter a slice of the Kerr geometry exists such that its three- 
geometry is conformally flat. It seems that, on the other hand, it is s impler to find solutions to the initial value 
problem compatible with Kerr perturbations on the usual hypersurfaces p3,p3 . In the present paper we will be more 
concerned with the question of how to impose initial data to the Teukolsky equation. 

Cauchy data in General Relativity consists of {gik,Kik) and a three-dimensional spacelike hypersurface Ef, where 
gik is a Riemannian three-metric on S^, and pl[ 



(2.3) 



is the extrinsic curvature which describes the "embedding" of Et in the spacetime, written in terms of the lapse 



~g'*)-i/2 and of the shift N, 



background metric gjlS , the Cauchy data for Einstein's equations, SG^!'} 



gti- In perturbation theory, to first order in the metric perturbation about the 



SnSTji],', are given in terms of the first 

order perturbations of the metric, , and in terms of the first order extrinsic curvature, on a constant time 

hypersurface. In the case of the Teukolsky equation, the first order initial value problem must be translated into 



terms of ^'|t=o and dt^\t=o- As we will show in this Section relating (/i|^', kH"^) at the initial hypersurface t 
(^t'|^=o, cJtVE'lt^o) is not straightforward. 

The starting point is given by Chrzanowski 



with 



I, who related the Teukolsky functions to metric perturbations 



-0^ ^ = ^US -3a- (3- Tr + f){5 - 2a - 2/3 + vr + f)/i„„ + {A-^i-fL-3j + j){A + fi-fl-2-/ + 27)/i„-,™ 
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2 {A-n-fl-3-f + 7)(f + 7T) + {5-3a- (3-Tr + f )(A - 2/2 - 27) /i,™ 



(2.4) 



and 



'^0^ ^ ^US + a + 3P - n + t){S + 2a + 2p - n - T)hii + (D + p + p + 3e - e){D + p - p + 2e - 2e)h„ 



2 {D + p + p + 3£-l){T + TT) + {5 + a + 3f3-Tr + T){D + p + 2e) h 



(2.5) 



where hnn = n^n'^hf^^, him — l^'m'^h^^jj^ etc, contain not only three-metric perturbations but all the components of 
the first order metric perturbations (here Z, n, m and rh are the null vectors of the Kinnersley tetrad). This represents 
a practical problem at the moment of imposing initial data, since one is normally given only the 3-geometry and the 
extrinsic curvature, not the 4-geometry and its time derivative (nor we should need it). The point is that ?/'4 (and 
"00 ) E^re gauge invariant quantities, but we do not have them explicitly expressed in terms of purely hypersurface 
quantities. 

Many simplifications in the analysis are possible when the background has spherical symmetry. In the Schwarzschild 
black hole case expressions (2.4)-(2.5) reduce to 



3 



4'' 



-^W [de - cote - 
16 \ r'^ \ smt 



dt ~ fdr + f 



+ ^ (dt - fdr + /') 



do- 



smt 



{htt - 2hrtf + Krf) 



2/ 



{dt - fdr) 



2i- 



T-9<p 



sm 



-cot9- 



smt 



Me 



smt 



- f [hre-i- 



sin 9 



(2.6) 



and 



(1) 



1 1 



de — cot 9 + 



sm t 



--(f-'dt + dr) id, 



cot ( 



sin( 



+ l^r'dt+dr + l^ {f-^dt+dr)^[l 



de + -^d^ {httf-^ - 2hrtf-^ + hrr) 



2i 



he. 



f-'{hte+i- 



sm 
h 



tip 



smt 



lire + i- 



(2.7) 



where / = 1 - 2M/r and /' = 2M/r'^ . 

The imposition of spherical symmetry carry also the following computational (but not necessarily crucial) advan- 
tage: the multipole decomposition of the metric perturbations in terms of spin- weighted harmonics -2Yl^{9) can be 
performed , and even and odd parity perturbations decouple so they can be considered independently. Below we 
shall decompose all metric perturbations in multipoles with (implicit) index £m (not to be confused with the tetrad 
vectors). 



A. Even or polar parity waves 



From now on we will consider only tp^ because of its better behavior for numerical integration compared to ipo 
and because the radiated energy can be described in terms of only ipA- Similar equations will hold for ipQ . Even 
parity perturbations are given by the real part (Re) of 1^4. In the Regge- Wheeler (RW) gauge, ho — hi = G = 



{hte = htif = hre = hnp 
derivative are given by 



0, hee = h^^/ sin 9), we have 



hw. 



0, then Eq. (|2.6D and its time 
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lor'' ^-^ 



(2.8) 
(2.9) 



Using Einstein's equations in the RW gauge, Zerilh's S, Eq. (C7e) gives. 
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(2.10) 



where i^fmis a source term given in Table III of Ref. jg]. The above expressions allow us to write Eqs. (2^) in terms 
of H2 and its time derivatives [Note that here and below there is an implicit index coming from the multipole 
decomposition of the metric perturbations.] The following step is to restore the gauge invariance, lost in principle 
when, for the sake of simplicity, we have chosen to work in the Regge- Wheeler gauge. To do so, we will use a result of 
Ref. where the metric perturbations in the RW gauge have been expressed in terms of 4>m , the Moncricf invariant 
function 
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which is expressed entirely in terms of the 3-geometry and is gauge invariant under changes of coordinates in the 
spacetime. The other piece of initial data, dt4>M , is expressed in terms of the extrinsic curvature Kij as in Ref. 
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where K^j^) and -^'(g) ^'^e respectively the "JC "and "G "parts of the extrinsic curvature, in analogy to the Regge- 
Wheeler decomposition of the metric tensor (see also a similar decomposition for K in Eqs. (2.32)-(2.35) of Ref. 
[^). Also, here and below, extrinsic curvature comp onent s (of even and odd parity) have indices tm implicit. 

In Ref. it was found by use of the definition (2.11) and the Hamiltonian constraint (sec for instance Zcrilli's 
equation (C7a) in Ref. ||^) 
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From the momentum constraint, equation (C7b) in Ref. |@ 
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The time derivatives of these metric coefficients in terms of (/) m have also been found in Ref. P6[ 
5tF2 = (r - 2M)dt(b'lj + , dtcf>'M 
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and from the time derivative of Eq. ( |2.14D , dtHi = rd^(j)'^,i + {XP - 3AMr - 3AP)/{r - 2M)/(Ar + 3M)d^(j)M , and 
the Zerilli [p| wave equation in vacuum, dfipM — d'^*(f>M ~ Vi{''')(f>M , we obtain 
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where Ve{r) is the Zerilli potential Here it is evident that in order to rewrite ipi and dtipA in terms of hypersurface 
data we had to use the evolution equations, not only the constraints. 

For the sake of simplicit y w e have written Eqs. (2.13)-( 2.16D in the sourceless case, but source terms can be 
straightforwardly included |2g]. The above relations together with ( 2.10 ) provide the expressions that allow us to 
write "04 and dtipi in terms only of 0m and dt4>M (and radial derivatives of them). This relation holds for all r and t 
and is gauge invariant since all involved final expressions are (in spite of having used the RW gauge as an intermediate 
step), in particular, this relation holds at the hypersurface where we want to give initial data. Since (pM can be wr itten 
only in terms of the perturbed 3-geometry and dt(t>M iu terms of only the extrinsic curvature (cfr. Eqs. (2.11) and 
( 2.12 )) , we will have, at St 



-04 = ■)p4:{hij,Kij), dtipi = dt')pi{hij,Kij) 



(2.17) 



our desired relations. Note that, in general, even for stationary, time- symmetric initial data neither tp^ nor dtip4 
vanish. This is in contrast with what happens with the corresponding initial data for the Moncrief wave form, where 
dt<f>M — 0, and with the initial data taken in Ref. |1^. 
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B. Odd or axial parity waves 



For odd parity perturbations we shall parallel the procedure we just followed in the even parity case. Let us recall 
that odd parity is described by the imaginary part (Im) of ip4 . In the RW gauge (/i2 — 0) we have 
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here hi and h^ are the odd parity metric perturbations (not to confuse with the same symbols used for the even parity 
perturbations, that in the RW gauge are taken to vanish). 

The corresponding gauge invariant waveform for odd parity given by Moncrief |9| is 
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Again, in the RW gauge (h^ = 0) we can write (from Eq. ( 2.20D ) 
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from Zerilli's B Eq. (C6c), one obtains 
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where D^rn. is a source term given in Table III of Ref . [pi , and finally /iq can be obtained from the definition of extrinsic 
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In spite of having constructed Eqs. ( p.22| )-( p.24| ) in the RW g auge , when we replace their right hand side in the 
expressions for Im'!/;4and lmdt'4'4, given by the Eqs. (2.18) and ( ^.19 ), they will be gauge invariant if we consider the 
general form of Q and K^p^ . We have then succeeded in our objective of expressing also the odd parity perturbations 
(the imaginary part of i/'4 ) in terms only of the 3-geometry and the extrinsic curvature. 



III. DISCUSSION 



In this paper we have faced (and solved for the a — Q case) the problem of giving initial data to the equation that 
describe in a gauge invariant way perturbations around the Kerr geometry, i.e. the Teukolsky equation in the time 
and frequency domains. We first observed that the connection between metric perturbations and gauge invariant 
objects that fully describe gravitational perturbations ipA (or V'o) as given by Chrzanowski |2^], is not explicitly 
written only in terms of the 3-geometry and the extrinsic curvature. This poses a difficulty when one wants to use 
this relation to impose initial data, since, normally one is given only the 3-geometry and the extrinsic curvature on a 
given hypersurface (and that should be all we need). In the nonrotating ( a = ) case (Schwarzschild background), 
we have been able to reexpress -04 (and V'o) and its time derivative in terms only of the 3-geometry and the extrinsic 
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curvature. To do so, we have momentarily gone to the Regge- Wheeler gauge and then restored gauge invariance by 
writing the metric perturbations in terms of (pM , the Moncrief's invariant waveform and its time derivative (and 
radial derivatives on the hypersurface) . As a by-product we obtained relations between the waveforms of Teukolsky 
and Moncrief (of course, they are not independent of each other) on any given hypersurface. These relations (for both 
the even and odd parity cases) represent the generalization of the Chandrasekhar transformations in the time domain 
when source terms are present. It is worth stressing here that to obtain the above described relations we used the 
evolution equations in order to eliminate the higher than the first time derivatives of the metric perturbations in the 
RW gauge. 

An interesting application of these relations is to compute the "close limit approximation" for black hole collisions 
by directly integrating the Teukolsky equation [|l9|. This constitutes a test for both, our relations on each hypersurface 
and the numerical accuracy of the code [ p^ for integrating the Teukolsky equation. 

All this could be accomplished for the nonrotating case, which provides the a — limit to be recovered by the 
more complicated (and realistic) case of a 7^ 0, for which many steps should be generalized. In fact, in the rotating 
case we do not have the multipole decomposition of the metric perturbations, nor a Regge- Wheeler gauge that can 
be reexpressed explicitly in terms of other invariants of the 3-geometry and the extrinsic curvature only, neither the 
generalization of for that case (that would be an interesting result in itself). For these reasons a more geometrical 
approach (making use of the Gauss-Codazzi relations) should be followed in order to rewrite ipi in the Kerr background 

A second aspect we have been able to study is the question of the regularization of the Green function solution 
to the Teukolsky equation in the frequency domain when nonvanishing initial data are present (see Appendix). The 
appropriate method to incorporate initial data in this case is to Laplace transform the Teukolsky equation. In this 
way, initial data appear as an additional source term and can be easily manipulated. We have explicitly shown in the 
nonrotating case and using Brill-Lindquist initial data that £ = 2 and £ = 3 modes need to be regularized. This has 
to be a general feature, valid also for the rotating background and for any other kind of astrophysically reasonable 
initial data for black hole collisions, and has to do with the fact that realistic initial data are not of compact support 
(or do not decay fast enough at spatial infinity). 
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APPENDIX A: REGULARIZATION OF THE INITIAL DATA 



In this appendix we shall consider again the problem of divergent integrals that appear when in the frequency 
domain one tries to solve the Teukolsky equation with sources that extend to infinity using the standard Green 
function technics. It was first proven by Poisson that, for a source generated by a particle released at rest from 
infinity (which has vanishing initial data), the divergent integrals appearing in the formal solution can be regularized 
in the nonrotating case. Soon after, the present authors |0 extended to the rotating (Kerr) background Poisson's 
approach. Here, we show how the same problem appears when one consider bounded sources, but nonvanishing initial 
data, like Misner or Brill-Lindquist data |29 3^ which are relevant to study black hole head-on collisions, in the close 
limit approximation. We thus apply the same regularization technics developed in references [|l6|,0 to show that the 
Green function formulation can be regularized. We shall refer to the above two references for further details of the 
problem and computations, and below only sketch how things go for the case of nonvanishing initial data. 

It is well known that the Teukolsky equation (2.1) can be separated by mode decomposition of the field and 
source into a complex radial wave equation and an angular equation satisfied by the so-called spin-weighted oblate 
spheroidal harmonics sSY^{0,aLL>) that in the case a = reduce ||32| to the spin- weighted spherical harmonics 
sY[^{0). By use of this separation of variables, the Teukolsky equation, in its frequency domain, was widely studied 
in the early seventies providing extraordinary results like the proof of the stability of rotating black hole and 

the computation of the gravitational radiation for unbounded particle trajectories around rotating black holes. 

A Fourier analysis in the complex plane can be used in the case one want to study the problem of a particle falling 
into the black hole from any finite distance. In this case one must consider the initial value problem, with ^> and 3f ^ 
specified on an initial hypersurface (at t = 0). To this purpose it is very useful to define the Laplace transform 'S' of 

to be 
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<b{uj,r,e,ip)= e''^'^{t,r,B,ip)dt . (Al) 

We take 4* to vanish for i < 0, which means that must be analytic in the upper half of the complex lo plane. 
Applying the Laplace transform to Eq. (2.1), the Teukolsky equation may be written in the following form 
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T(w, r, 6,ip)= / e*'^*47rST(t, r, 6*, 93) dt 
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In this way Cauchy data are explicitly incorporated into the Teukolsky equation as part of the source term. Eq. ( [A2| ) 
can be separated into a radial and an angular part pO| , p^ and applied to study particular problems involving finite 
radii infall. 



Let us consider the source term given by Eq. (A3) in the nonrotating case 

1 



SlD 
A2" 



r4 1 



2A/^ 



[{iujr^ + 4r - 12M) iP^it = 0, r) - r'^dtM* = 0, r)] ^Gg 



(A5) 



where G = (I - + + 2)(2^ + l)/(47r) , and g{r) takes the following form for the Brill-Lindquist initial 
data (30) 



4r - lOA/) 



(r - 2My 



1/2 



(A6) 



where r is the coordinate in the conformal space, that we can take as the isotropic radial coordinate, related to r as 
~= (\/^ + ^r-2hlf /4 . 

This expression holds for the case of two equal masses black holes initially at rest, and in the perturbative regime 
(close hmit [|,||), but a very similar expression holds for the case of unequal masses (included the particle limit). In 
particular, the conclusions about regularization will hold the same also for different initial data such as the Misner 
ones 1^ (symmetrization through the throats). 

As noted by Poisson |^ the divergent integral generated by the Green method, /div , can be regularized by 
introducing a function h(r) and integrating by parts i. e. writing it in the following form 



Gg{r)gi^e'^'CX 



dr 



H,oo\ 



dr - he^^'CX 



iut n \rH.OO 



(A7) 



boundaries 



where ^div = '^iujr'^ (1 - 3M/r + icor) , C = d/dr* + iuj (with r* = r + 2M\n {r/2M - 1) ), and X"^°^ are solutions 
to the Regge- Wheeler equation [ pO[ (in the frequency domain) with purely ingoing and purely outgoing behavior in 
the horizon and at infinity respectively. The function h[r) is determined by 



8 



h = -Ge- 



gGdivdr. 



By direct power counting in the first addend of Eq. (A7) we observe tliat the more important modes for computing 
gravitational radiation, modes £ = 2 and ^ = 3 (in the unequal masses case) need to be regularized since the integral 
diverges linearly and logarithmically, respectively. 



Proceedings of the November 1994 meeting of the Grand Challenge Alliance to study black hole collisions may be obtained 

by contacting E. Seidel at NCSA (unpublished). 
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